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Abstract 

The maximum magnetic efficiency of a generai muitipoie Haibach cylinder of order p is found as function of p. 
The efficiency is shown to decrease for increasing absoiute vaiue of p. The optimai ratio between the inner and 
outer radius, i.e. the ratio resuiting in the most efficient design, is aiso found as function of p and is shown to tend 
towards smailer and smailer magnet sizes. Finaily, the demagnetizing field in a generai p-Halbach cylinder is 
caicuiated, and it is shown that demagnetization is iargest either at cos2/70 = 1 or cos2/70 = -1. For the common 
case of a p = 1 Haibach cyiinder the maximum vaiues of the demagnetizing fieid is either at ^ = 0, tt at the outer 
radius, where the fieid is aiways equai to the remanence, or at 0 = ±n/2 at the inner radius, where it is the 
magnitude of the fieid in the bore. Thus to avoid demagnetization the coercivity of the magnets must be iarger 
than these vaiues. 
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1. Introduction 

Generating a powerful magnetic field with a permanent mag¬ 
net assembly is important in a number of applications. It 
is also often of importance that the magnetic field is gener¬ 
ated using the least amount of magnet material possible. A 
common magnet design used to generate such powerful mag¬ 
netic fields is the Haibach cylinder, which has been used in a 
number of applications, such as nuclear magnetic resonance 
(NMR) equipment [1; 2], accelerator magnets [3; 4], mag¬ 
netic refrigeration devices [5; 6] and medical applications 
[7]. The Haibach cylinder is a hollow cylinder made of a 
ferromagnetic material with a remanent flux density which in 
cylindrical coordinates is given by 

^rem.r — ^rem COS ptj) 

Sin p(j), (1) 

where Z?rem is the magnitude of the remanent flux density and 
p is an integer [8; 9]. For p positive a field inside the cylinder 
bore is generated, which for the important case of p = 1 is 
spatially uniform. There is zero field outside the cylinder. For 
p negative the Haibach cylinder creates a field on its outside, 
while the inside field becomes zero. Usually, as will also 
be the case here, a two dimensional problem is considered. 
This is a good approximation if the radius of the cylinder is 
smaller that its length. The magnetic field distribution for such 
a Haibach cylinder of infinite length [10; 11; 12; 13; 14] as 
well as for finite length [15; 16; 17; 18] have previously been 
investigated in detail. However, the efficiency of Haibach 
cylinders have not been considered for a general p-Halbach 
cylinder. 

The plan of the paper is as follows: First we calculate the 
efficiency of a p-Halbach cylinder and thereafter we consider 



Figure 1. The Haibach cylinder geometry for a p-Halbach 
ip > 0). The different radii and regions, I-III, have been 
indicated. 


the demagnetization field internally in the Haibach cylinder 
and its possible influence on the performance and efficiency 
of the Haibach cylinder. Finally, we discuss the implications 
of our findings. 

We consider the field from a p-Halbach in general. The 
geometry is as shown in Fig. 1 for the case of an ‘interior’ 
Haibach (positive p); for an ‘exterior’ Halbachs (correspond¬ 
ing to negative p) the field is generated on the outside of the 
cylinder. We assume that the permanent magnets are perfectly 
linear and with a relative permeability of =1, which is a 
good approximation for NdFeB magnets. 
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2. Efficiency of a Haibach cyiinder 


The object of a permanent magnet array is to generate a mag¬ 
netic field of given characteristics in a given volume. Many 
different magnet configurations can in principle produce the 
same magnetic field, and thus the question arises of how to do 
it most efficiently. Jensen and Abele [19] proposed a general 
figure of merit, M* , to characterize the efficiency of a given 
magnet design; 




where y^eid is the volume of the region where the magnetic 
field is created and Vmag is the volume of the magnets. The 
figure of merit is the ratio of the energy stored in the field 
region to the maximum amount of magnetic energy available 
in the magnetic material. It can be shown that the maximum 
value ofM* is 0.25 [19]. 

We here consider the efficiency of a p-Halbach cylinder. 
While the efficiency of a p = 1 Haibach cylinder have been 
considered several times [18; 20; 21], the efficiency of a 
general p I Haibach cylinder have not been considered in 
the literature, although the calculation is straight-forward. The 
efficiency can be found using the analytical solution to the 
field equations for the magnetic field for a general p-Halbach 
cylinder. These are given in Bjprk et. al. [14], where the field 
interior to the cylinder is given as: 



R,/R 

I 0 

Figure 2. The efficiency, M*, as function of the ratio of the 
inner and outer radius for Haibach cylinders with p = —5 to 
p^5. 


It is of interest to determine the optimal efficiency possible 
for a general Haibach cylinder and the ratio of the inner and 
outer radius at which this occurs. This can be done by taking 
the derivative of M* (Eq. (5)) with respect to the ratio of the 
radii and equating it to zero. The resulting equation is, as can 
be seen below, a polynomial equation of order p+\, which 
does not in general have a closed form solution for p > 3: 
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Here Rg is the outer radius of the Haibach and Ri the inner 
radius, as shown in Fig. 1 . 

For an exterior Haibach (p < —1) the field outside the 
Haibach is; 
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Performing the integral in Eq. (2) for an internal and ex¬ 
ternal Haibach cylinder yields 
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The efficiency as a function of the ratio between the radii 
is plotted in Fig. 2 for p = —5 to p = 5. 


For the case of p = 1 the solution is given by ^ ^)/2- 

where W is the Fambert W function. The argument of the 
Lambert W function is greater than — 1 /e, which means that 
the function is single-valued. Evaluated numerically this cor¬ 
responds to the well known ratio of ^ = 2.2185, although 
only the numerical solution has previously been given in lit¬ 
erature [20; 21]. For other values of p Eq. (6) can be solved 
numerically. The solution, i.e. the highest efficiency and cor¬ 
responding ratio of the radii, is shown in Figs. 3 and 4 as 
a function of p. In order to verify the analytical results for 
the field, the value of M* for the optimal ratio of the radii 
was also calculated numerically using Comsol Multiphysics 
finite element software (FEM) and was found to match the 
analytical value. 

For the case of p = —1, the Haibach cylinder is simply a 
uniformly magnetized circle or infinite rod. For this system 
the efficiency becomes: 



This clearly has a maximum for Rj = 0, which means there is 
no internal bore. In this configuration the efficiency is 0.25, 
which as noted is the maximum efficiency possible. 
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Figure 3. The optimal ratio between the inner and outer 
radius, i.e. the ratio producing the highest efficiency M*, as 
function of p. 
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Figure 4. The highest value of the efficiency for a p Haibach 
cylinder as function of p. 


3. Demagnetization effects 


When a given Haibach design is implemented using perma¬ 
nent magnets with a finite coercivity, the demagnetization 
field internally in the magnet becomes an issue that need to 
be considered. Inside the magnet, the magnetic field H is 
not in the same direction as the magnetization M, and if the 
component of the magnetic field in the direction opposite to 
the magnetization exceeds the coercivity. He, of the magnetic 
material, the magnet will be demagnetized [22], with severe 
implications for the produced magnetic field and therefore also 
for the magnetic efficiency of the magnet design. Therefore it 
is important to consider demagnetization effects together with 
the magnetic efficiency. For the p = \ Haibach it is known 
that the reverse component of the magnetic field can exceed 
the intrinsic coercivity in regions around the inner equator 
[17; 26], but for the general p Haibach the demagnetization 
field is not known. 

Using the expressions for the field inside the magnet from 
Ref [14] it is straightforward to calculate the demagnetization 
of a Haibach cylinder. The condition for demagnetization to 


occur is 

H M 

- <-Hc. (8) 

^rem 

We consider the cases p = 1 and p > 1 separately. The case of 
p < 1 is not considered, but the calculations follow the case 
of p > 1. 


3.1 Demagnetization in a p = 1 Haibach 

The vector field inside the Haibach cylinder is given by [14] 

= Beem{rinRo - r\nr)sm^. (9) 


From this equation the magnetic flux density is readily de¬ 
rived in the same way as in Ref. [14], and using elementary 
trigonometric relations we get 

=M(r^Biem ( (In^ - ^)cos20 - ^ ) , (10) 
£)rem \ ^ J 


The condition for demagnetization becomes 


^rem 


T^Xln 



r 


-)cos20. 


( 11 ) 


Values of (r, (j )) which fulfil the above inequality corresponds 
to the places in the Haibach cylinder where demagnetization 
will occur. As long as the prefactor on the cosine is pos¬ 
itive, demagnetization will first arise for cos20 = — 1, i.e. 
0 = ±;r/2. This is on the ‘equatorial plane’ of the Haibach 
cylinder. If the left hand side of Eq. (11) is positive, i.e. if 
poHc/Brem < 1 /2, then for any value of the inner radius R, 
there will be demagnetization. 

Let us then consider poHc/Brem > 1 /2. Demagnetization 
will always occur first at r = R,. As R,- is decreased there will 
be a critical value for which demagnetization first occurs on 
the inside of the Haibach cylinder. This critical value is given 
by 


PoHc Ro 

-= In- 

^rem ^m.crit 


( 12 ) 


or 


R 


m,crit — exp 


^rem / 


(13) 


This is exactly equal to the value of the internal field produced 
by the Haibach cylinder. Demagnetization will thus first occur 
where the magnetization is pointing in the opposite direction 
of the field in the bore, and the value of the demagnetization 
field will be the Haibach field due to field continuity across 
the border between magnet and bore. 

If the prefactor ln(Ro/r) — 1/2 is negative, demagnetiza¬ 
tion will first occur for cos 20 = 1, corresponding to 0 = 0 or 
7Z. In this case, demagnetization will occur from the critical 
radius: 


RoHc Ro 

-= 1 — In- 

^rem ^crit 


( 14 ) 
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Figure 5. A polar plot of the demagnetization field for a 
p—I Haibach cylinder with Ro = l,Ri — 0.2 and Brem = 1 T 
for a coercivity of /ioHc = 0.8 T. 


or 



to Ro, with the largest value occurring at r = Ro, where the 
demagnetizing held will always have a value of Brem- 

Thus, for 1/2 < jj^Hc/Brem < 1 demagnetization can oc¬ 
cur hrst both at 0 = ±;r/2andat(/) =0,7:, while for/roHc/Brem 
1 demagnetization always occur at (j) = ±;r/2 hrst. In other 
words, the demagnetizing held is always Brem at 0 = 0, ;r and 
r = Ro while it is BremIn ^ at ^ = ±7r/2 and r = B, and these 
are the maximum values of the demagnetizing held and thus 
to avoid demagnetization the coercivity must be larger than 
these values. An illustration of the demagnetization held for a 
p = l Haibach cylinder with Ro = I, Ri = 0.2 and Brem = 1 T 
for a coercivity of p^Hc = 0.8 T is shown in Fig. 5. 

3.2 Demagnetization in a B>1 Haibach cyiinder 

The calculation of the demagnetization for p > I proceeds 
very similarly to the case of p = 1. Using the equations of 
the internal held in the magnet from Ref. [14] we get that the 
condition for demagnetization to occur becomes 


BqHc 1 / P f 

Brem 2 P - I \Ro 




lB+l\ 

2p-lj 


coslptj). 


( 16 ) 


When the prefactor of the cosine is positive (i.e. when r < 
Bo(l/2 + l/2p)'/(P^*)), demagnetization arises hrst for 
cos2p0 = — 1. Again, there will always be demagnetiza¬ 
tion for poHc/Bi-gm < 1/2 no matter the value of B, . For 
PoHc/Brem > 1 /2, the critical inner radius is given by 


Only for materials with an He small enough that the right 
hand side of Eq. (17) is positive, can demagnetization occur. 
It is evident that higher p cylinders are less susceptible to 
demagnetization. 

For r>Bo(l/2-|-l/2p)*/(f’^') demagnetization hrst arises 
for cos2p0 = 1, and in this case the demagnetization from Bo 
to the critical radius is 

1,P-1BoHo 

This can only occur if I /2 < poHc /B^m < 1 ■ 

Similarly to the p = 1 case, for p > 1 there will always 
be critical demagnetization inside the Haibach cylinder for 
poHc/Brem <1/2. For 1/2 < BoHc/Bi-em < 1 the hrst demag¬ 
netization can occur both for cos2p0 = 1 and cos2p0 = —1, 
while for poHe/Beem > 1 demagnetization always occur hrst 
at cos2p(/) = —1. If BoHe/Brem > P/(p ~ 1) there will be no 
demagnetization at all in a p Haibach cylinder. 

The demagnetization held calculated above is for Haibach 
cylinders with a continuously rotating remanence. In most 
practical applications a Haibach cylinder is segmented, which 
can cause demagnetization to occur also at e.g. the interfaces 
between segments [17]. However, the demagnetization held 
as calculated above will still generally be present and must 
be taken into account also for segmented designs. Several 
Haibach cylinder-like high held systems (> 3T) have been 
constructed using permanent magnets [23; 24; 25]. Here de- 
> magnetization is to some extend avoided by replacing magnets 
subject to a high demagnetization held with a high permeabil¬ 
ity soft magnetic material. However, the constructed magnets 
still produce a lower Hux density than expected, which can be 
caused by the demagnetization issues discussed above. 

4. Discussion and conciusion 

Here we have calculated the magnetic efficiency of a general 
p-Halbach cylinder and showed that the efficiency decreases 
for increasing absolute value of p. We have also found the 
optimal ratio between the inner and outer radius, i.e. for 
the most efficiency design, as function of p. The optimal 
ratio cannot be expressed in closed analytical form, but the 
numerical solution is given. Finally, we have consider the 
demagnetization effect in a general p-Halbach cylinder, and 
showed that demagnetization is largest either at cos2p0 = 1 
or cos2p0 = — 1. For a p = 1 Haibach cylinder the maximum 
values of the demagnetizing field are at 0 = 0, ;r and r = 
Ro where the field is always Brem, while it is the magnitude 
of the field in the bore at (j) = ±7r/2 and r = B,. Thus to 
avoid demagnetization the coercivity must be larger than these 
values or the magnets at these locations must be replaced by a 
high permeability soft magnetic material. 
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